Wind and structural loads data measured on parabolic trough solar collectors at an operational power plant

Wind loading is a primary contributor to structural design costs of concentrating solar-thermal power collectors, such as heliostats and parabolic troughs. These structures must resist the mechanical forces generated by turbulent wind, while the reflector surfaces must maintain optimal optical performance. Studying wind-driven loads at a full-scale, operational concentrating solar-thermal power plant provides insights into the wind impact on the solar collector field beyond the capabilities of wind tunnel tests or state-of-the-art simulations. We conducted comprehensive field measurements of the atmospheric turbulent wind conditions and the resulting structural wind loads on parabolic troughs at the Nevada Solar One plant over a two-year period. The measurement setup included meteorological masts and structural load sensors on four trough rows. Additionally, a lidar scanned the horizontal plane above the trough field. In this study, we describe the high-resolution dataset characterizing the complex flow field and resulting structural loads. This first-of-its-kind dataset will enhance the understanding of wind loading on collector structures and will help in designing the next-generation solar collectors and photovoltaic trackers.

deflections can significantly impact the optical performance of solar collectors, thereby reducing the overall plant efficiency.
To generate a comprehensive dataset of wind loading on collector structures, we performed an extensive measurement campaign at an operational parabolic trough power plant.The aim was a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs to provide insights into structural dynamic response.NREL initiated the field measurement campaign at the Nevada Solar One (NSO) parabolic trough plant (https://solarpaces.nrel.gov/project/nevada-solar-one).Operating since 2007, the plant has a nominal capacity of 72 MW with 0.5 h storage, and uses thermal oil as the receiver working fluid at temperatures up to 393 °C.Constituting the solar collector field, parabolic troughs are aligned in the north-south direction and track the sun from east to west over the course of a day, facing varying wind loads depending on the wind properties and the angle of the troughs.The collectors (type SGNX-1) are designed as an aluminum space frame construction with glass mirrors.During strong winds and nighttime, the troughs are stowed 30° below the horizon to the east.
The measurements at NSO, running from October 2021 through June 2023, were recorded by sonic anemometers on masts at different heights to characterize the incoming flow and conditions around four trough rows at the western edge of the trough field.In addition, a Doppler lidar scanned the horizontal plane above the troughs.These wind measurements were combined with structural load measurements.The load measurements were installed on the same four outermost trough rows and included support structure bending moments, drive torque moments, dynamic accelerations of the space frame, mirror displacement, and tilt angles.Our measurements stand out for several reasons.First, they encompass a combination of wind conditions across multiple trough rows, along with corresponding loads and the wind field directly above the trough field.Second, our measurements extend over a long time period, allowing us to capture a broad variation in environmental conditions.Last, we employed high-frequency measurements, enabling us to study the dynamic interactions and intricate dynamics of the system in detail.The present paper describes the extensive field measurements of the wind masts, structural loads, and lidar.In addition to using these data for designing next-generation solar collectors, the dataset can be used to create and validate computational models for predicting the unsteady flow conditions and wind loading in collector arrays.

Methods
Overview of the measurements at Nevada Solar One.The NSO parabolic trough plant is located near Boulder City, Nevada, USA, at 35.8 • N, −114.983• E and at 540 m elevation in a hilly desert environment.The plant's solar collector field consists of 95 loops, each with eight 100 m long solar collector assemblies (SCAs) aligned in the north-south direction.NREL's measurement setup is located at the western edge of the collector field at the four outermost rows, targeting winds perpendicular to the trough rows.Figure 1 shows an overview of the instrumentation setup.Meteorological masts for in situ wind measurements and load sensors on the collector structures are strategically installed within 49 m long half-sections of an SCA (Fig. 2).The trough angle of each collector assembly is controlled by a drive actuator in the center of the assembly (the northern end of our measured segment), situated on the drive occurrence (DO) pylon.The loose bearing at the southern end is referred to as the shared occurrence (SO) pylon.More details about the trough geometry are given in the section about structural loads measurements.The lidar is located 70 m south of the inflow mast outside of the field.Figure 3 summarizes the time periods of data availability for the different instruments.The masts continuously operated between October 2021 and June 2023.Due to regular maintenance on the parabolic troughs between April 2022 and November 2022, the NREL team decommissioned the wake towers between the rows in this time period.Structural loads measurements were added in November 2022.The lidar operated throughout the entire period in different scanning configurations, but only the periods of published data are shown in Fig. 3.
Mast measurements of winds.General setup.The mast setup consists of an inflow mast outside of the array and three wake masts between the four westernmost trough rows.Table 1 and Fig. 3 provide an overview of the individual instrument characteristics on all masts and their operating period.The inflow mast with sonic anemometers at heights of 7 m, 5 m, and 3.5 m characterizes the incoming flow on the western side of the trough field.The sonic anemometers provide wind speed components in the north, east, and upward directions as well as sonic temperature with a temporal resolution of 20 Hz.The 7 m sonic anemometer measurements extend  above the parabolic troughs, and the 3.5 m sonic anemometer is located approximately at the hinge height of the troughs.A cup anemometer adds the horizontal wind speed at 15 m height.At 2 m height, temperature, relative humidity and barometric pressure were measured at 1 Hz frequency.The configuration was slightly changed over the course of the measurement period: The 5 m sonic anemometer was added first in February 2022, the lowest sonic anemometer had to be replaced twice due to availability of the instruments, and the temperature sensors at 3.5 m and 7 m were first added in November 2022.
Additionally, three wake masts between the four westernmost rows characterize the wind field between the trough rows with sonic anemometers at heights of 3.5 m, 5 m, and 7 m.Again, the configuration was slightly modified due to practical reasons-most importantly, the height of the lowest sonic anemometers was changed from 4 m to 3.5 m in February 2022.The 5 m sonic anemometer on mast 3 was not installed in the period starting November 2022.
Prior to installation, all meteorological instruments were calibrated.The cup anemometer and Gill WindMaster sonic anemometers were calibrated in an accredited wind tunnel, while the CSAT3 sonic anemometers and Gill WindMaster HS sonic anemometers (new units) were calibrated at the manufacturer.The NREL calibration lab, an ISO 17025-accredited facility, calibrated the Vaisala PTU (pressure, temperature, and humidity) sensors in an environmental chamber and the Rosemount temperature sensors in an oil bath.
A Campbell CR1000X data logger collected the digital input signals of all met mast instruments and a GPS receiver.The GPS time stamp was used for synchronization with other data streams.Data were transferred to NREL servers hourly using the data logger server software via a cellular modem.The program's primary scan, at 20 Hz, was used to poll the sonic anemometers with secondary scans at 1 Hz for the GPS synchronization, cup anemometer, temperature, humidity and pressure sensor.The GPS ensured the real-time clock accuracy was within ± 10 μ s.Raw data were sent to an NREL server and stored as hourly files separately for the inflow 20 Hz and 1 Hz data and for each wake mast.These files were then used for postprocessing.
Data processing.In the data postprocessing routines, the hourly raw data files are read and merged into daily datasets for the inflow mast and combined wake masts.The general processing steps for all sonic anemometers are identical.For all parameters, outliers are excluded by detecting values outside a 5-times standard deviation (5σ) in a 60 s rolling window.The amount of detected outliers is consistently below 2% of daily data for all sonic anemometers.The wind coordinate system in the raw sonic data, which is based on the specific sonic model, is transformed into a right-handed local tangent plane coordinate system with wind speed components u E in the east direction perpendicular to the parabolic trough rows, v N in the north direction along the trough rows, and upward w ("ENU").This system is also consistent with the loads coordinate system.All Gill sonic anemometers we affected by the "w-boost" firmware bug 17 , which includes an underestimation of measured vertical wind speeds.We apply the correction factors of 1.166 for upward velocities and 1.289 for downward velocities to the raw data.Then, the horizontal wind speed U is calculated by and the meteorological wind direction WD (indicates where wind is coming from) results from with respect to the quadrant of u E and v N in the defined coordinate system.The sonics in the field were oriented with a calibrated digital protractor; no further tilt corrections were made due to the complexity of the surrounding environment.The low-resolution meteorological data, including temperature T, pressure p, and relative humidity (RH) at 2 m height, are added to the sonic wind speed data.The sonic output includes high-frequency measurements of sonic temperature (which is approximately equal to virtual temperature); however, these sonic temperature values are not calibrated.No consistent calibration has been found for the time period when actual temperature data at the sonic height levels were available.Despite this limitation, the sonic temperature measurements still capture high-frequency variations in temperature, which are used for covariance estimations (see below).The separate temperature sensors installed at the inflow mast yield more accurate absolute temperature values, although with a larger response time.Additionally, for each sonic anemometer, the turbulence parameters turbulence kinetic energy (TKE) and turbulence intensity (TI) are calculated in a running 10 min time window for each time step, ensuring a high resolution of TI and TKE.The TKE is determined by summing the squares of the velocity standard deviations in the three orthogonal directions: where σ uE , σ vN , and σ w represent the standard deviations of the wind velocity components.The TI along the mean horizontal wind U, for the vertical wind component w, and for the wind components perpendicular to the troughs u E and along the troughs v N , is calculated as the ratio of the standard deviation of the respective wind speed component to the mean wind speed: , TI , TI and TI (4) These turbulence parameters provide insights into the dynamic characteristics and intensity of turbulence within the measured wind field and in relation to the array of parabolic troughs.
At the inflow mast, additional parameters are derived to characterize the state of the atmospheric boundary layer.To assess atmospheric stability, the bulk Richardson number is calculated between the 3.5 m and 7 m height level, for time periods with temperature measurements at these levels.The bulk Richardson Number (Ri b ) is defined by: where g = 9.81 m s −2 is the acceleration due to gravity, T ref is the average temperature between the two levels, and Δz = 3.5 m is the height difference with differences in temperature ΔT and in the horizontal wind speed components Δu E and Δv N .As part of the surface energy budget, turbulent energy fluxes describe the energy transport by turbulent eddies within the atmospheric boundary layer.These fluxes are directly determined from our high-frequency measurements using covariances in consecutive 20 min time intervals.The (virtual) sensible heat flux H S , derived from the covariance of (virtual) temperature and vertical wind speed, quantifies the turbulent heat transfer between the surface and the atmosphere.The momentum flux τ, derived from the covariance of horizontal and vertical wind speed, characterizes the transfer of momentum due to wind shear.These turbulent fluxes are computed as: The overline denotes an average of the subrecord, while the prime signifies the turbulent fluctuation ) .The mean air density ρ is calculated based on the 2 m inflow measurements, and c p = 1006 J kg −1 K −1 is the specific heat capacity of air.The covariance w U ′ ′ includes fluctuations from all horizontal wind com- ponents as suggested by Weber 18 for the friction velocity.Covariances are computed for parameters at the 7 m height level from detrended, non-overlapping 20 min time segments.The 7 m height level was selected since this measurement height was least impacted by the solar collector field (the maximum collector height is ∼5.5 m), thereby providing the most representative measurements of undisturbed atmospheric wind conditions.Ogive analyses 19,20 have demonstrated that a 20 min window effectively captures the major frequency contributions to the covariances and confirm their convergence within this time frame for the majority of cases.Covariances also allow for the calculation of the Obukhov length and the flux Richardson number, which are alternative measures of dynamic stability based on high-frequency measurements at only one height.The Obukhov length L is defined by the vertical turbulent fluxes of heat and momentum u U w is the friction velocity as calculated for the momentum flux and κ = 0.4 is the von Kármán constant.The 2 m temperature serves as T ref , since this variable is available for the entire measurement period.The sign of L indicates the stability of the atmosphere: positive values correspond to stable conditions and negative values to unstable conditions.The flux Richardson number R f is determined by with the covariance of the eastward wind component u E and vertical wind u w E ′ ′ and the wind gradient u z / E Δ Δ between 3.5 m and 7 m (analog for v N ).For all sonic anemometers on the wake and inflow masts, turbulent integral length scales are also calculated.Integral length scales L describe the characteristic size of energy-containing eddies in a turbulent flow.Here, they are calculated from the autocorrelation function A of a continuous time series from a point measurement (as provided by the sonic wind speeds).A w for the vertical wind component w(t) results from where A(t * ) represents the autocorrelation function at time lag t * .The integral of the autocorrelation curve first provides the integral time scale T w , which can be transferred into a length scale by L T = ⋅ U w w using Taylor's hypothesis of "frozen turbulence" transported by the mean wind 21 .The detrended 20 min time series of U and w are shifted by time lags t * between 0 s and 20 min, and for each lag, the autocorrelation coefficient A of the original with the shifted time series is calculated.As an approximation for the integral under the A w (t * ) curve, the curve is assumed to be an exponential function, and the intercept yields the integral time scale 19,21 .If the time scale is larger than the time series segment, the result is excluded.We find that the exponential fit method yields the most robust results, although estimating length scales from time series is challenging and many large length scales are missing because of the limited length of the time series segments.The same approach as discussed for w is applied to derive the turbulent length scales L U for the horizontal wind speed and uE L and L vN for the eastward and northward component, respectively.The postprocessed measured and derived data are stored at 20 Hz resolution with the covariance-derived values added at 20 min intervals.Variables in the processed dataset are listed in Table 2. Additionally, 1 min averaged data are provided with peak values and standard deviations of wind speeds in these intervals.
Lidar measurements of winds.To characterize flow conditions directly above the solar collector field, a Galion G4000 pulsed Doppler wind lidar 22 was installed approximately 70 m south of the inflow meteorological (met) tower.This scanning lidar was installed on top of two shipping containers at approximately 6 m above the ground (Fig. 1) with a leveled top surface.We used a Ushikata Transit to align the Lidar in the field along the cardinal compass directions.
The Galion lidar works by emitting laser pulses into the atmosphere and analyzing the backscattered light to measure the Doppler shift, providing information about line-of-sight (LOS) wind speed in the direction of the laser beam.Along the beam, the lidar provides a spatial resolution corresponding to the non-overlapping, 18 m range gate lengths.Usable data can be obtained between ≈ 70 m distance from the lidar (the "blind zone") to the configured maximum scanning distance of 612 m.The lidar at NSO was configured to do two types of scans in the horizontal plane above the collectors: a 360 ° horizontal scan (or Plan Position Indicator-PPI scan, used between June 2022 and September 2022) and a 0.4 Hz stare scan (or fixed LOS scan, between April 2023 and May 2023).The purpose of the 360 ° horizontal scan was to investigate the impact of parabolic troughs on wind conditions within the plant at a height above the troughs.This pattern performed full horizontal 360 ° scans in 1 ° increments every 10 minutes and yielded reliable data within an opening angle of ± 30 ° about the mean wind direction.During a stare scan, the lidar operates and scans at a fixed azimuth and elevation angle, and the fixed laser beam captures temporal variations at a higher temporal resolution.The stare scan in the west-east direction over the solar collector field aimed to characterize turbulence conditions along the west-east scanning direction during western wind events.Using the ≈ 0.4 Hz LOS wind speed within each range gate, the variance in wind speed can be determined.Note that resulting variances might be underestimated due to the limited probe volume averaging.Further, western wind directions are rare at NSO, limiting available data for the stare scan period, despite scientific interest.The technical aspects of the two different scan patterns are listed in Table 3.
The dataset (variables are listed in Table 4) includes the LOS wind speed along the laser beam ("Doppler") and information about the direction of the beam, characterized by azimuth (Az, the angle in the horizontal plane starting at north) and elevation (El, the angle above the horizontal plane).For the horizontal scan, Az varied between 0° and 360° and El was constant at 0°, whereas both Az and El were constant at 90° and 0°, respectively, for the stare scan.Before analyzing wind speed data, low-quality data need to be filtered by calculating the signal-to-noise ratio (SNR) and eliminating data samples with low SNR.The SNR is calculated from the power intensity PI (provided in the lidar data) using: For our data, we recommend using an empirical threshold of SNR = −19 such that lidar data with SNR > −19 is kept for further analysis while the rest is discarded.However, the threshold might have to be adjusted for specific cases.The filtered LOS wind speed is included in the dataset ("Doppler_filtered").Note that pitch and roll represent the lidar instrument alignment.These values were not changed during the course of the measurement campaign, and no alignment correction was made in the dataset.Lastly, the lidar dataset contains some short interruptions leading to multi-hour gaps in data during the summer months due to the lidar overheating.
Measurements of structural loads.General setup.In addition to the meteorological mast measurements, structural loads at the parabolic trough collectors were measured between November 2022 and June 2023.The sensors were installed in rows 1, 2 and 4 on the half-SCA close to the met masts.The individual load sensors and their processing routines are explained in more detail in the following sections.The sensor locations on the space frame or trough mirrors are depicted in Fig. 4. The individual sensor models are listed in Table 5.In addition to loads sensors, a GPS time stamp was sampled as well as the cup anemometer wind speed signal at 15 m height.The data from all load sensors were collected in one central data acquisition system (DAS).The DAS is based on highly configurable National Instruments (NI) hardware and software.A local host computer in the plant control room and a NI PXI system (composed of chassis, controller, and modules) were used to control and interface with three chassis that were distributed across trough rows 1, 2, and 4. Communication between the PXI system and the chassis was established through fiber optics, employing EtherCAT media converters.This setup allowed for reliable and efficient data transfer between the components.As a first signal processing step, calibration slopes and offsets were applied to the raw data in the DAS so that the resulting files contained data in engineering units, with the exception of the strain gage signals where the conversion factors were applied in postprocessing.The collected data were saved to the DAS PC.Data were regularly uploaded from the PC to a file transfer protocol site hosted by NREL and then routinely downloaded and saved to the NREL server for postprocessing.During the test campaign, there were several periods where loads data were lost, probably due to interrupted communication from the PXI to the EtherCAT data chain.These events required the DAS to be manually reinitialized to resume communications and data collection.The data acquisition process involved a data scan rate at 1 kHz, while data were sampled and stored at a reduced rate of 20 Hz, applying an anti-aliasing filter in the NI signal conditioning.
For postprocessing, the loads data had to be time-synchronized with the measured wind data.This was done with the GPS time stamp recorded in both data streams.For loads measurements before February 23, 2023, the loads GPS sensor did not receive satellite signals.Therefore, the time stamp of the data logger had to be corrected manually by correlating the wind speed signal at 15 m height, which was recorded in both the inflow mast and loads datasets.This manual time stamp correction provides an accuracy of below 1 s for the affected period.As a first general data processing step, outliers and nonphysical values were excluded.Variables in the processed dataset are listed in Table 6.All published loads data are provided as time series at 20 Hz resolution and additionally at an averaged 1 min resolution.These sampling frequencies match the met mast frequencies and enable straightforward merging of the datasets.Additional parameters in the 1 min loads dataset are maxima, minima, and standard deviations for all loads variables for users interested in peak values or dynamics.The original 20 Hz time series can also be used to calculate loads spectra.Individual sensors.Support structure bending momentsBending moments of the trough support structure were measured at rows 1, 2, and 4 at both the SO and DO pylons.The bending moments were measured using strain gages in a full bridge T-configuration that is shown in Fig. 4a.The strain gages measure deformations that are calibrated to provide an experienced moment.The DO and SO pylons are built differently (Fig. 2): the DO pylon has a triangular truss construction and the SO pylon is a simple bending beam; therefore, the installation of the gages differed at the DO and SO pylons.For both, the gages were placed near the lower third of the exposed pylon and at least one width of the tube/web away from the concrete foundation.For the drive pylon,  (c) Inclinometer near shared occurrence.
Fig. 4 Pictures of the load sensor installations.
the gages were placed at H gage = 0.480 m from the foundation and on the inside of the c-channel.This location was selected to avoid interference from electrical enclosures attached to the outside of the pylon.The shared occurrence gages were placed at H gage = 0.175 m from the foundation and on the inside of the square tube.At this location, the gages were approximately one cross section away from the base plate.
To arrive at engineering units of the strain signals, the strain signal offset was subtracted from the raw strain data.The strain offsets were determined using the average strain values over the entire dataset, filtered for low wind speeds and bin-averaged over tilt angles less than ±70° to the upward position.Next, the signal was multiplied by the scale factor.The scale factors were determined through a calibration pull of the instrumented structural component, using a telehandler with slings, shackles, and a reference load cell in line with the rigging (Fig. 5).Load was incrementally applied by slowly retracting the boom.This calibration provided moment values for each SO and DO pylon.
To convert the bending moments at the measurement location to overturning bending moments M B (the actual bending moment about the foundation axis parallel to the y direction), the following correction was applied with H c = 0.279 m (Fig. 2): The bending moments can also provide a drag force in the x-direction: The measured moments are defined right-handed positive along the y-axis (Fig. 2).Because of the different layouts of the pylons, bending moments behaved differently at SO and DO.At the simple-bending-beam SO pylon, the measured moments can be assumed to be actual bending moments.In contrast, at the triangle-truss-structure DO pylon, the measured moments show some overlaying torque moment characteristics (resulting from tension and compression in the truss structure) in addition to pure bending characteristics.As a result, the bending moments and drag forces at the DO should be treated with caution, but are still included in the dataset for completeness.
To compare the forces and moments to other studies, the loads were transformed to load coefficients C mb for the bending moments and C fx for the drag forces: where the 60 s mean air density ρ and mean horizontal wind speed U were determined from the met mast data- set.L segment = 8 m is the spanwise length of a mirror segment, and = W 5 m is the aperture width.The moment coefficient C mb represents the bending moment about the foundation of the pylons, whereas the horizontal drag force coefficient C fx is calculated from the measured bending moment by dividing by the lever length (hinge height).The load coefficient definitions are adapted from Hosoya et al. 9 .Coefficients were calculated for data with wind speeds (inflow mast at 3.5 m height) greater than 3 m/s to avoid excessive values resulting from low values in the denominator.For 1 min averaged data, the coefficients were calculated with the 1 min means of moments, wind speed, and air density.The load coefficients for wind speeds from eastern directions should be treated with caution because the wind field is impacted by the entire trough field before reaching the instrumented rows, and the coefficients might not be reasonable.For the sake of completeness, they are still included in the published data.Drive torque momentsTorque moments at the drive location were measured on rows 1, 2, and 4 using strain gages configured as a full torque bridge.The strain gages were placed on the driveshaft south of the drive actuator (Fig. 4b).There was limited free length of the drive shaft to allow for strain gage installation.The gages were instead installed on the shaft of the hub that attaches to the drive shaft.In this case the installation did not meet the recommended distance of one cross-section diameter of the drive shaft away from any possible stress riser or discontinuities.The gages were oriented to 90° away from the shaft keyway to avoid discontinuities.Torque moments were calibrated in the field analog to the bending moments using the mean of low-wind, tilt-bin-averaged offsets and a calibration pull for slope calibration.
The torque moments are defined right-hand positive about the pivot axis y of the parabolic troughs (Fig. 2).The torque moment distribution over the torque axle has its maximum at the drives (corresponding to the measured moment) and decreases linearly to zero near the SO end of the torque axle.Similar to the bending moment, the torque moments have been transformed into moment coefficients C my : where L panel = 49 m is the spanwise length of an entire mirror panel (one-half of an SCA).Dynamic tiltDynamic tilt angles were measured at row 1, row 2, and row 4. Each instrumented trough row included three dynamic inclinometers located near the bottom of the space frame with one unit at the SO (Fig. 4c), one unit at the DO, and one centrally located ("Mid") on the trough section.The inclinometers were installed so that when the trough is flat, the response is a zero-degree reading.Positive degrees are measured when the trough is facing east (in the morning) and negative degrees when facing west (in the evening).The stow position is 120 ° (30 ° below the eastern horizon), but occasionally the troughs are brought into a stow position of 150 ° or 90 °.The trough's "nominal" tilt angle can be calculated from the sun angle according to the NREL Sun Position Algorithm (https://midcdmz.nrel.gov/solpos/spa.html),which is available through a python API 23 .The nominal tilt angle θ nom is calculated through the algorithm described in Anderson & Mikofski 24 , which projects the sun position into the transversal (east-west) plane using the elevation and azimuth angles.The dataset includes this transversal sun angle β EW at the NSO location.
Measured tilt angles θ 0 at horizontal trough positions for each sensor were calibrated in the field using digital protractors.However, the inclinometers measure the angle of the space frame at the measurement location instead of the actual angle of the parabolic troughs.Comparing measured tilt angles against expected trough angles (from the sun angle), we found consistent offsets.We provide a corrected tilt angle that accounts for these offsets.For this, we define the tilt angle error ε along the transversal (east-west) plane as ε = θ 0 − θ nom .Using a known angular position-stow at 120°-we calculate an average calibration constant for each sensor in the 1 min data and subtract this constant from the data.To find this calibration constant, we extract data that satisfy the following conditions: • Measured tilt angle equals 120° ± 1°.
• Wind speeds are less than 1 m/s.• Maximum tilt angle error ( 120 ) is less than 1.6° because tilt angle errors greater than this value indicate operational settings outside of stow.
We apply this filtering to the data and calculate average tilt angle errors for the time periods before and after December 22, 2022, separately because the NSO operators performed a calibration on this date and recalibrated the drive motors of the parabolic troughs.
The filtered dataset consists of 5,159 data points before December 22, 2022, and 3,747 data points after December 22, 2022.The average tilt angle errors are calculated for each sensor and are shown in Table 7.We subtract the average tilt angle error for each sensor from the tilt data before and after December 22, 2022, to calculate the calibrated tilt angle θ θ ε θ 0 n om , which reduces the average tilt angle error.For example, in row 1 at the DO, the average tilt angle error of all daytime data after December 22, 2022, and below the error threshold of 1.6° is reduced from −1.51° to 0° after the calibration.The dataset contains both "raw" tilt angles as measured by the inclinometers θ 0 as well as calibrated values θ with the average stow position error subtracted (Table 6).Mirror deflectionsMirror deflections were measured using laser displacement sensors on two single mirrors on the western edge of trough rows 1 and 4. Each of the two instrumented mirrors is located at the center span of the trough section south of the intermediate pylon.Five laser displacement sensors at each mirror were placed on the space frame behind the mirror to measure the distance between the rear surface of a mirror near each of its corners (NW, SW, NE, SE; in proximity to the connection points) and at the mirror center (Fig. 4d).The sensors attached to the space frame provided the relative distance between the mirrors and the space frame.As a result, the measured displacements represented a combination of the support structure deformation and the deformation of the glass mirror itself.This relative distance was between 75 mm and 80 mm for the central sensors and between 98 mm and 108 mm for the sensors at mirror corners (Fig. 6).To find a "zero-position" for the displacement measurements, we filter cases with low winds (<3 m/s) and vertically oriented mirrors (to exclude gravity influences as shown in Fig. 6a).This vertical mirror position occured in the late mornings before solar noon.Unfortunately, the opposite position (180° rotated) was never taken in normal operation.We noticed that the described zero-position is not a constant displacement value but depends on the ambient temperature.The relation of relative displacement to temperature is nearly linear (Fig. 6b) with a positive slope for the corner sensors and a slightly negative slope for the center sensors, but with a spread below 1 mm for all sensors.We define a zero-displacement at the vertical mirror position and at an ambient temperature of 15 °C to study wind, temperature and gravity influences on mirror deformations.The dataset contains both "raw" relative distances between mirror and space frame as well as adjusted values with the zero-displacement subtracted (Table 6).Accelerations at the space frameAccelerometers were installed on rows 1, 2, and 4 on the northern and southern end of one-half SCA near the DO and SO pylons.They were attached to the space frame on the western edge of the trough (Fig. 4e).The accelerometer sensors were attached to the spaceframe behind the mirrors and rotate with the assembly.The biaxial measurement directions of each accelerometer are orthogonal (with the local y-axis in the plane of the mirror and with the local x-axis perpendicular to the mirrors).The frequency range allowed for a wide range of dynamics measurement, down to the static acceleration, or gravity.The static acceleration due to gravity changes from 1 g to 0 g as the sensor moves from a vertical orientation (aligned with gravity) to a horizontal orientation, and therefore also provides a measure of the trough angle as this is proportional to the tilt.The acceleration signals for low-wind cases were corrected to be within +1 g and −1 g, when the trough was rotating.These corrections were subsequently applied to all data irrespective of the wind speeds.the column structure of the Parquet format allows to selectively read only the desired columns.Example Python scripts to read the datasets in Parquet file format are provided at the OEDI repository 25 .The data structure consists of daily Parquet files, where separate files are provided for the 20 Hz and 1 min resolutions (for loads and met masts).The naming convention for these files follows the pattern: Type_resolution_YYYY-MM-DD_00h_ to_YYYY-MM-DD_00h.parquet,where YYYY-MM-DD represents the date of the daily file that contains data from 00:00 UTC to 00:00 UTC the next day (times may vary when there are data gaps).The Type can be one of four options: Inflow_Mast, Wake_Masts, Loads, or Lidar, while the resolution can be either 1 min or 20 Hz (the lidar files contain no resolution).All time stamps in the data are in UTC format.For example, the file name for 1 min loads data on November 19, 2022, would be Loads_1min_2022- 11-19_0h_to_2022-11-20_0h.parquet.For the lidar data, each file includes one complete lidar scan, and the file names include the hour, minute, and second of start and end time.An example is Lidar_2023-04-02_19-00-17_to_2023-04-02_19-00-17.parquet.For the lidar stare scan, only the data for 15 selected western wind periods are archived, whereas continuous lidar data are available for the horizontal scan period (June 2022-September 2022).Table 2, Table 4, and Table 6 list all variables in the met mast, loads, and lidar datasets with their respective formula symbol, a description, and units.When a single instrument on a met mast was not available on a specific day, the respective data column is filled with NaNs.In addition to these listed variables, the 1 min met mast datasets include maximum and standard deviation values of wind speed, derived in 1 min windows.Furthermore, the loads datasets contain maximum, minimum, and standard deviation values of 1 min windows for all loads variables.

technical Validation
Mast measurements of winds.The met mast measurements at NSO are compared to the METAR (METeorological Aerodrome Report) observations from Boulder City Municipal Airport, located approximately 20 km northeast of NSO.These data are provided by the Iowa Environmental Mesonet of Iowa State University (https://mesonet.agron.iastate.edu/request/download.phtml).Note that differences in observations can arise due to spatial, orographic, and temporal variations.However, this comparison serves to verify the accuracy of the coordinate transformations applied to the sonic anemometers.Figure 7 presents a comparison of wind speed and direction measurements across all met mast heights.Additionally, the temperature and RH measurements from the inflow mast are included in the comparison.The results indicate a good qualitative agreement between the met mast measurements at NSO and the METAR observations.This agreement lends confidence to the correctness of the coordinate transformations used for the sonic anemometers and validates the reliability of the measured wind speed, wind direction, temperature, and RH data.
Lidar measurements of winds.The wind speeds determined by the lidar are verified with the inflow met tower measurements.Since the lidar measures wind speed at about 6 m above the ground, the measurement at 7 m height from the inflow met tower is used for qualitative validation.The LOS wind speed from the lidar is first converted into true wind speed using the wind direction from the inflow met mast before conducting the comparison.Overall, the lidar-interpreted wind speed matches qualitatively well with the measurements from the inflow met mast (coefficient of determination r 2 = 0.85; Fig. 8), suggesting good agreement between these two Fig. 7 Verification of the met mast measurements at Nevada Solar One with METAR data from the Boulder City Municipal Airport.independent instruments.Note that r 2 in this case is not greater than 0.9, which is a common threshold lidar-mast comparisons.This is most likely because the distance between lidar and inflow met tower is too close to the measurement blind zone of the lidar (four range gate numbers, or approximately 72 m).More lidar-associated analysis will be presented in a follow-on study.

Usage Notes
An example Jupyter notebook script written in python is provided at the OEDI repository 25 to read the presented datasets.

Fig. 1
Fig. 1 Overview of wind and structural loads measurements at NSO. Satellite images: © 2023 Google Earth Data.

Fig. 2
Fig. 2 Bending and torque moment definitions at the drive occurrence (DO) and shared occurrence (SO) with trough dimensions and the loads/wind coordinate system.The sketch shows one half of an SCA.Image courtesy Besiki Kazaishvili, NREL.

Fig. 3
Fig. 3 Data availability throughout the measurement period at NSO.

( d )
Laser displacement sensors visible on the NW, SW and center locaƟon of the mirror.(e) Accelerometer at western edge of mirror.(b) Torque gage placement at the drive shaŌ.Half of the torque bridge is shown.(a) Drive pylon bending gage locaƟon and installed T-bridge gages (leŌ) and shared post T-configuraƟon strain gage installaƟon (right).

Fig. 5
Fig. 5 Schematic and pictures of strain gage calibration pulls.

Table 1 .
Instrumentation for the inflow and wake masts.

Table 2 .
Variables included in the published inflow and wake mast datasets.

Table 5 .
Variables included in the published lidar datasets.Instrumentation for the structural loads measurements.

Table 3 .
Technical details of the two types of lidar scan patterns.

Table 6 .
Variables included in the published loads datasets.